Non-monotonic temperature dependence of thermopower in strongly correlated electron systems 
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We examine the temperature dependence of thermopower in the single band Hubbard model using dynamical- 
mean-field theory. The strong Coulomb interaction brings about the coherent-to-incoherent crossover as tem- 
perature increases. As a result, the thermopower exhibits non-monotonic temperature dependence and asymp- 
totically approaches values given by the Mott-Heikes formula. In the light of our theoretical result, we discuss 
the thermopower in some transition metal oxides. The magnetic field dependence of the thermopower is also 
discussed. 
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Thermopower is none other than the amount of entropy flow 
along with the electric current. 1 The consideration on the en- 
tropy in thermodynamics tells us the low and high temperature 
(T) limits of the thermopower: In the metallic systems, the 
thermopower goes to zero as T — > 0. On the other hand, the 
high-temperature limit of thermopower is given by the entropy 
consideration 13 in the atomic limit. In the strongly correlated 
systems, the spin and orbital degrees of freedom enhance the 
high-temperature thermopower^^ 

In the low-temperature limit, the ratio of the thermopower 
Q and T is proportional to the derivative of the density of 
states D{uS) with respect to the energy uj measured from the 
chemical potential p as Q/T oc —dD(uj)/dui\ ul —o. By this 
relation, not only the sign but also the magnitude of the ther- 
mopower of conventional metals is well understood. The ther- 
mopower is a sensitive tool for the electronic states. 

The electron correlation brings about exotic electronic 
phases such as an anomalous metal near the Mott transition. 
In the vanadium oxide, Lai-^Sr^VOs, the filling control Mott 
transition is realized, and non-monotonic temperature de- 
pendence of the thermopower is observed.. 4 . The temperature 
dependence manifests a crossover of coherent-to-incoherent 
charge transport. This phenomenon is common to transition 
metal oxides. The cobalt oxide, Na x Co02 is an example^— 
In the photoemission spectroscopy measurementSjiSiii it is 
reported that the coherent motion of charge carriers is 
rapidly suppressed with increasing temperature, and the quasi- 
particle-peak disappears at ^200 K. 

In this paper, we study the role of the strong Coulomb in- 
teraction on thermopower, whose temperature dependence is 
particularly examined in detail. For this purpose, the single 
band Hubbard model is adopted as a minimum model and the 
strong Coulomb interaction is treated in the dynamical mean 
field theory (DMFT) 12 , which can capture the coherent-to- 
incoherent crossover due to the strong Coulomb interaction. 
This method based on the local picture is useful to understand 
the overall behavior of thermopower as a function of temper- 
ature. We find that the Coulomb interaction significantly af- 
fects the temperature and magnetic field dependence of the 
thermopower. The Coulomb interaction is found to give rise 
to a non-monotonic temperature-dependence, which is well 
described by the entropy consideration at high temperatures. 



In the light of our theoretical results, we discuss the thermo- 
electric response in some transition metal oxides. 
We start with the single band Hubbard model: 
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where er(=t,|) denotes electron spin, is the dispersion 
relation of the non-interacting electrons and other notations 
are standard. In the DMFT formalism, the resulting equa- 
tions are the functions of the density of states for the non- 
interacting electrons. We denote the "bare" density of states 
by Do, and take the semicircular function of the energy e, 
D {e) = [2/(ttW 2 )]VW 2 - e 2 , with the normalization W = 
1 throughout this paper. Because the Hamiltonian has the 
particle-hole symmetry, all the results shown here are in the 
case that the electron concentration n > 1. To solve the sin- 
gle impurity problem in DMFT, we employ the non-crossing 
approximation (NCA)riii and the iterated perturbation the- 
ory (WT) i i — which do not require the analytic continua- 
tion from the imaginary frequency axis. Using the calculated 
spectral density p ff (e, u) through DMFT, the thermopower Q 
is expressed by, Q = — {ks / e){Ai/ Aq), where 
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where p a {e,ui) = Im[l/{w + p — e — S CT (w)}] with 
the electron self-energy. 

In the earlier studies, a saturation behavior of ther- 
mopower at high temperatures in the Hubbard model was 
discussed^— We find in the following that the asymptotic 
behavior of Q at high-temperatures shows a non-monotonic 
temperature dependence. Before presenting the numerical re- 
sults, let us note the high-temperature thermopower Q of the 
model Eq. (Q~|i, on the entropy consideration. The indepen- 
dent variables of the function Q are the electron concentration 
n, the Coulomb interaction U, and temperature T. For fixed 
n, we have two high-temperature limits: i) Q\ := Q(T — >• 
oo, U) by keeping fcgT < U. Because U —> oo is achieved 
before T -> oo in this case, Q\ is given by, 
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ii) Qi ■= Q{T -> oo, U) by keeping k B T > t/. In this limit, 
U is of less importance, so that the result is written as the well 
known Heikes formula, 
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Therefore, we expect the two different asymptotic behaviors, 
i.e., Qi and Q 2 , and furthermore, a sign-change of the ther- 
mopower may occur in the temperature dependence: Figure 
[TJa) shows the n dependence of the high-temperature limits 
of the thermopower, Q\ and Qi. For n < 1.3, Q\ is negative 
whereas Q 2 is positive. 

Figure [Tib) shows the temperature dependence of the ther- 
mopower calculated by DMFT with the NCA impurity solver. 
We find the non-monotonic temperature dependence of the 
thermopower. This is well understood as the asymptotic be- 
havior with the high-temperature limits, Qi and Q 2 : In the 
temperature region, fc^T > 0.4, Q is in the range between Qi 
and Q 2 for each n (see the closed and open dots in Fig.QJa)). 
With increasing T, Q approaches Qi first, and further increas- 
ing T, Q shows the saturation behavior given by Q 2 . For 
n = 1.4, Q is always positive. But for n = 1.3, Q changes its 
sign twice with increasing temperature. With further decrease 
in n, the absolute value of the minimum of Q is enhanced and 
becomes closer to Q\. This is because the Coulomb inter- 
action U is more effective near half filling. The effect of the 
Coulomb interaction is made much clear by the U dependence 
in Q. In Fig.QIb), the results of U = 4 and 7 are shown for the 
same electron concentration, n = 1.1. We see that the asymp- 
totic approach of Q to Qi is obvious for larger U. It is worth 
noting that the entropy consideration on the high-temperature 
thermopower works well even at finite temperatures. 

The single band Hubbard model Eq. (fl]i has 2-fold degen- 
eracy on the singly-occupied site. This is a disadvantageous 
condition for the NCA impurity solver because the approxi- 
mation is based on an expansion in 1/N where N is the ionic- 
angular-momentum degeneracy: For ksT < 0.4, Q increases 
rapidly with decreasing temperature. The numerical calcula- 
tion simultaneously becomes unstable and eventually breaks 
down at certain temperature, i.e., the imaginary part of the 
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FIG. 1: (Color online) (a) Thermopower in the high temperature 
limit, Qi and Q2, vs. electron density, n. (b) Temperature depen- 
dence of the thermopower calculated by DMFT with the NCA impu- 
rity solver, for various parameter sets of U and n. 




FIG. 2: (Color online) Temperature dependence of the thermopower 
for n = 1.2 calculated by DMFT with the IPT impurity solver. Bro- 
ken and dotted lines are Qi and Q2 for n — 1.2, respectively. 



calculated self-energy becomes positive at small frequencies. 
This will be improved for large N systems, though further 
considerations are required for the thermodynamic properties 
near zero temperature 

The thermopower is sensitive to the Coulomb interaction 
even for small U. To examine the thermopower in the small U 
region, we use the IPT impurity solver for DMFT. In the cal- 
culation, we follow the method by Refs. Il8lfl9ll to obtain the 
self-energy. Figure |2]shows the temperature dependence of Q 
for n = 1.2 and various U. In the small U limit, U = +0, 
we use Dq(cj + fx) for J dep a (e,uj)D (e) in Eq. ©. In this 
case, Q is a monotonically increasing function of ksT: In the 
low temperature limit, Q/T oc — dDo(ui + fj,) / dLo\ u —o and at 
high temperatures, Q asymptotically approaches Q 2 , i.e., the 
Heikes formula Eq. (@}. On the other hand, Q shows the non- 
monotonic temperature dependence for finite U. Near zero 
temperature, the gradient of Q with respect to T is larger for 
larger U . With increasing T, Q shows a maximum and then a 
minimum, and increases with further increasing T. For large 
U, the convergence of the numerical calculation in DMFT 
with the IPT impurity solver becomes poor, and the minimum 
of Q grows beyond Q\ (see the result for U = 3.0 in Fig. 2). 

We note the complementary characters of the impurity 
solvers NCA and IPT in DMFT. The care must be taken to 
discuss the results of the thermopower shown in Figs. [T]and 
[2] The NCA impurity solver is a perturbative expansion in 
powers of the hybridization between the impurity site and the 
effective bath. Therefore, NCA is appropriate to discuss the 
thermopower for large U (we find that the results are not plau- 
sible for U < 3). On the other hand, the IPT impurity solver 
is a perturbative expansion with respect to the Coulomb in- 
teraction U and then appropriate for small U. Further as an 
advantage, IPT can access the very low temperatures unlike 
NCA. Next, we discuss the response of the low-temperature 
thermopower to the magnetic field by DMFT with the IPT 
impurity solver focusing on the relatively small Coulomb in- 
teraction U < 2W (total band width of the non-interacting 
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the total density of states, D tot (w), is written as D{oj + B) 
D(lu — B), the derivative dD tot (uj)/duj is expressed as 
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FIG. 3: Magnetic field dependence of the thermopower for fesT = 
0.02 and n — 1.2 calculated by DMFT with the IPT impurity solver. 



system). 

We introduce the Zeeman term, — BJ2i( n it ~ n ii)< mto 
the Hamiltonian Eq. (Q~|i to discuss the response of the ther- 
mopower to the magnetic field B. Figure [3] shows the mag- 
netic field dependence of the thermopower for kgT = 0.02, 
n = 1.2 and several values of Coulomb interaction U. We 
find that the Coulomb interaction enhances the thermopower 
and its response to magnetic field. The result for U = +0 is 
also presented as a reference to discuss the role of U on the 
response of Q to B. For finite U, the convergence of the nu- 
merical calculation becomes poor. Therefore, in Fig. [3] we 
plot Q within the range of B where IPT is accessible. 

In general, Q in the low temperature region is well ex- 
plained by the structure of density of states around the Fermi 
level. Hence, the increase of Q by B is also understood from 
the same viewpoint (see Fig. @). In the small magnetic field, 
we expand the interacting density of states D a (oj, B) with re- 
spect to B as, 

Du(w, B) =D{u±B)=Y. ^r-^(± B ) m - (5) 



m=0 



duo r 



Note that D t (u,B = 0) = D±(uj,B = 0) = D(lu). Since 
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FIG. 4: (Color online) The density of states (DOS) for n = 1.2 at 
k B T = 0.02. Thick line is the result of DMFT with the IPT impurity 
solver for U = 1.5, 5 = 0. For reference, DOS for U = 0, B = 
is presented by a thin line. 



dD tot {u:) 
dw 
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Through the relation Q/T oc — dDtot{u) / du\ u= o for the 
low-temperature thermopower, the magnetic-field dependence 
reflects the detailed structure of the density of states D(oj) 
near u> = 0. For [7 = 0, the condition, semicircular density 
of state D with n = 1.2, results in an increase of Q for small 
magnetic field, and the slow response to B is a consequence 
of the non-interacting electron system. A simple dome struc- 
ture of Do manifests the increase in Q against B. Figure [4] 
shows that the density of states for U = 1.5 has a coherence 
peak near ui = 0. Therefore, as seen in Fig. [3] the increase of 
Q also appears. 

Equation © suggests that the response of the low- 
temperature thermopower to magnetic field reflects the de- 
tailed structure of the density of states near the Fermi level, 
i.e., the increase or decrease in the thermopower under the 
magnetic field is dependent on the differential coefficients, 
d 2m+1 D tot (u)/d(jj 2m+1 \ u =o. In this theoretical study, we 
use the single band Hubbard model with semicircular density 
of states Do. The thermopower within DMFT is a function of 
the density of states. Therefore, the increase in Q by the mag- 
netic field as shown in Fig. [3] is a consequence of the simple 
model where density of states has a negative slope at u = 
for n = 1.2. In reality, however, reflecting the band structure, 
the density of states near the Fermi level is certainly differ- 
ent from the simple dome shape. Here, the more important 
is that the response to the magnetic field is enhanced by the 
Coulomb interaction U. As shown in Fig.|4]by the thick line, 
the coherence peak is created by the Coulomb interaction ac- 
companying the abrupt change in the density of states near the 
Fermi level. As a result, the response of the thermopower to 
the magnetic field is enhanced as shown in Fig. [4] 

In the present study, the qualitative behavior of the ther- 
mopower is clarified in the wide range of temperature, al- 
though we use the simplest model. An essential feature 
demonstrated here is that the strong Coulomb interaction 
brings about the large response of thermopower to external 
disturbance through the instability of the electronic state char- 
acterized by the narrow coherence peak. A measure of the in- 
stability is the renormalized energy scale, i.e., the width of the 
coherence peak. The renormalized or reduced energy scale is 
of importance not only for the magnetic field dependence but 
also for the temperature dependence of the thermopower. In 
the temperature dependence, the entropy consideration for the 
high-temperature thermopower Q\ works even at much lower 
temperatures than the band-width 2W and/or the Coulomb- 
interaction U (see Figs.QJa) and[2]). This means that the high- 
temperature is achieved on the basis of the reduced energy 
scale by the Coulomb interaction. 

In the vanadium oxide, (La,Sr)VC<3, a non-monotonic tem- 
perature dependence of the thermopower is observed^ The 
thermopower approaches the two limiting values expected 
from the entropy consideration. This is the evidence of the 
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coherent-to-incoherent crossover of the electronic states and 
well explained by our theory. In the cobalt oxide, Na^CoC^, 
a strongly renormalized quasi-particle band, which disappears 
near the room temperature, is observed by the photoemis- 
sion spectroscopy measurements As discussed in this pa- 
per, when electrons lose their coherence by increasing tem- 
perature, the thermopower simultaneously approaches to the 
asymptotic value obtained by the entropy consideration in 
the correlation dominant regime, i.e., Qi for fc^T < U. 
Hence, for the large thermopower of this material observed 
near the room temperature, the strong Coulomb interaction 
must be one of important factors. For more qualitative stud- 
ies, one should employ the developed DMFT analysis such 
as cluster-DMFT 24 and (cluster-)DMFT combined with the 
density functional theory.— In fact, the recent cellular DMFT 
studies with realistic band structures suggest that the non-local 
correlations improve the results.— Yet, we believe that our 
study provides a good starting point for the future theoretical 
studies based on the advanced methods. On the other hand, in 
the partially-occupied t2 g states, the spin-orbit coupling is not 
fully quenched as discussed in Refs.l27landl28l In electronic 
systems with heavier elements such as S^IrOi, the spin-orbit 
coupling becomes comparable to the Coulomb interaction and 
the kinetic energy and, thus, plays a significant role for the 



Mottness of correlated systemsj22r— This effect may be an- 
other path to a new thermo-electronics based on the strongly 
correlated electron systems. 

In summary, we have theoretically studied the role of the 
Coulomb interaction in the thermopower. To clarify the role, 
we consider the single-band Hubbard model, and calculate 
the thermopower using the dynamical mean field theory. We 
find the non-monotonic temperature dependence of the ther- 
mopower. This is well described by the asymptotic form in 
the high-temperature limit given by the entropy consideration. 
The Coulomb interaction plays an important role to obtain 
such asymptotic behavior in a rather low-temperature region 
by creating the narrow quasi-particle band. 

After completing the manuscript, we have noticed a recent 
work by W. Xu et al£ who studied the high-frequency limit 
of the thermopower in the strongly correlated system. They 
also reported the non-trivial temperature dependence. 
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